r  fiO-A181  193  THE  FRAM  CAVITATION  CODEOJ)  UEIDLIMGER  ASSOCIATES  NEU  1/1 
VORK  N  REHfiK  ET  AL  15  HAV  86  UA-86-7  DNA-TR-86-179 
0NA881-84-C-8881 

UNCLASSIFIED  F/G  12/5  ML 


ioaaaMXiowMtta^nonrww^,yrou«RX7VXiui  ?u  rh  m  vxlt  mnwawwjBwgwiwanumuiin  ot 


AD-A181  193 


DNATR-86-179 


THE  FRAM  CAVITATION  CODE 


M.  Rehak 
R.  Smilowitz 
R.  Kagel 

Weidlinger  Associates 
Consulting  Engineers 
333  Seventh  Avenue 
New  York,  NY  10001 


15  May  1986 


Technical  Report 


CONTRACT  No.  DNA  001-84-C-0001 


Approved  tor  public  release; 
distribution  is  unlimited. 


THIS  WORK  WAS  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODE  B344085466  Y99QMXSF00045  H2590D. 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  DC  20305-1000 


87  6  8  ^6 


t'j  *•* 


6  6J 


VV  yjy  '.J_  JTLjn^KUHjjrM.  jvww~-w*  i--h  i.h  *r»t  xtutm 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return 
to  sender. 

PLEASE  NOTIFY  THE  DEFENSE  NUCLEAR  AGENCY 
ATTN:  TITL,  WASHINGTON,  DC  20305  1000,  IF  YOUR 
ADDRESS  IS  INCORRECT,  IF  YOU  WISH  IT  DELETED 
FROM  THE  DISTRIBUTION  LIST,  OR  IF  THE  ADDRESSEE 
IS  NO  LONGER  EMPLOYED  BY  YOUR  ORGANIZATION. 


CUT  HERE  AND  RETURN 


DISTRIBUTION  LIST  UPDATE 


This  mailer  is  provided  to  enable  DNA  to  maintain  current  distribution  lists  for  reports.  We  would 
appreciate  your  providing  the  requested  information. 

□  Add  the  individual  listed  to  your  distribution  list. 

□  Delete  the  cited  organization/individual. 

□  Change  of  address. 

NAME:  _ 

ORGANIZATION:  _ 

OLD  ADDRESS  CURRENT  ADDRESS 


TELEPHONE  NUMBER:  _i _ I_ 

SUBJECT  AREA(s)  OF  INTEREST: 


DNA  OR  OTHER  GOVERNMENT  CONTRACT  NUMBER: _ 

CERTIFICATION  OF  NEED-TO-KNOW  BY  GOVERNMENT  SPONSOR  (if  other  than  DNA): 
SPONSORING  ORGANIZATION:  _ 

CONTRACTING  OFFICER  OR  REPRESENTATIVE: _ 


Director 

Defense  Nuclear  Agency 
ATTN:  Hi  TITL 
Washington,  DC  20305-1000 


Director 

Defense  Nuclear  Agency 
ATTN:  M  TITL 
Washington,  DC  20305-1000 


UNCLASSIFIED 


t 


SsC.R'ty  C.A$S.=  CA~ON  of  °ags 


dMiim 


REPORT  DOCUMENTATION  PAGE 


a  R£FOR’  SECURITY  CLASS.F  CATION 
UNCLASSIFIED 


•O  RESTRICTIVE  MARK.NGS 


2a  SECURITY  CLASSIFICATION  AUTHORITY 

N/A  since  Unclassified 


Lb  DECLASSiFiCATlON ,  DOWNGRADING  SCHEDULE 
N/A  since  Unclassified 


3  DISTRIBUTION  i  AVAILABILITY  OF  REFORT 

Approved  for  public  release; 
distribution  is  unlimited. 


4  FERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

WA  86-7 


5  MONITORING  ORGANIZATION  REPORT  NUMBEP(S) 

DNA-TR-86-179 


6a  NAME  OF  FERFORMING  ORGANIZATION 

Weidlinger  Associates 
Consulting  Engineers 


6b  OFFICE  SYMBOL 
(If  applicable) 


7a  NAME  OF  MONITORING  ORGANiZAT.ON 

Director 

Defense  Nuclear  Agency 


6c.  AOORESS  City,  Start,  ana  ZIP  Coat) 

333  Seventh  Avenue 
New  York,  NY  10001 


7b  ADDRESS  (City,  Surf,  ana  ZIP  Code) 
Washington,  DC  20305-1000 


3a.  NAME  OF  finding  ,  SPONSORING 
ORGANIZATION 


8b  0FF.CE  SYMBOL 

(if  applicable) 
SPSS/Tsai 


3  FROCLREMEN*  INSTRUMENT  DENTiFiCAT'ON  NUMBER 

DNA  001-84-C-0001 


3c  AOORESS  (C/ty,  State  ana  ZIP  Coat) 


'0  SOURCE  OF  F'jNQiNG  NUMBERS 


PPOG^AW 

PRO. EC* 

~A$< 

WORK  _N;T 

Element  no 

NO 

NO 

ACCESS  ON  NO 

62715H 

Y99QMXS 

F 

DH008838 

' '  riTL£  (include  Stcunfy  Classification) 
THE  FRAM  CAVITATION  CODE 


•2  PERSONAL  AuTHORtS) 

Rehak,  M. ;  Smilowitz,  R.  and  Kagel,  R. 


’3a  ’*voE  OP 

13b  time  COVERED 

•4  DATE  OF  RERCR’  rear  Month  Oay) 

15  face  COUNT 

Technical 

=rom  851031  -o 8605 15 

860515 

32 

6  supplementary  notation 

This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMSS  Code  B344085466 
Y99QMXSF00045  H2590D. 


COSATI  COOES 

’8  SuBjECT  terms  (Continue  on  reverse 

f  necessary  ina  identify  by  block  number ) 

=  ElD 

GROUP 

SUB-GROUP 

•Bulk  Cavitation 

-.Underwater  Explosions, 

13 

LO 

Characteristic  Analysis 

19 _ 

_ 2 _ 

Surface  Ship  Response^ 

'9  ABSTRACT  Continue  on  reverse  •  f  necessary  ana  'dentify  by  blocs  number I 


This  report  describes  the  computer  code  FRAM.  The  code  computes  the  pressure  and 
velocity  in  a  fluid-half  space  where  bulk  cavitation  may  occur.  The  cavitation  results 
from  an  exponentially  decaying  pressure  wave  from  an  underwater  explosion  which  travels 
to  the  surface  and  is  reflected  as  a  tension  wave.  Because  the  water  does  not  sustain 
tension,  a  cavitated  region  is  formed.  The  cavity  subsequently  closes  under  gravity 
(lower  closing  shock)  and  accretion  of  vapor  particles  (upper  closing  shock). 

The  theory  is  summarized  herein.  Capabilities  and  limitations  are  discussed  and 
a  description  of  the  source  code  is  included. 


20  DISTRIBUTION'  AVAILABILITY  OF  abstract 
□  UNCLASSiF  ED'UNl.MiTEQ  CT  save  as  RpT 
22*  NAME  OF  RESPONSIBLE  NOiviOUAl 
Sandra  E.  Young 

OD  FORM  1473,  34  mar 


□  PTC 


jSERS 


21  A8STRAC'  SECLRi’y  C.ASSiF  CAT  ON 

UNCLASSIFIED 


22b  telephone  (include  Area  Code)  22c  OF=CE  svN1BOl 

(202)  325-7042  DNA/CSTI 


33  APR  taition  may  oe  ..lea  jiti  jirjutr ea 
AH  art  QOSO'tt# 


CuASS^'CA^QN  Qp  ^H,s  SAG 
UNCLASSIFIED 


CONVERSION  TABLE 


Conversion  factors  for  U.S.  Customary  to  metric  (SI)  units  of  measurement. 


MULTIPLY 
TO  GET  — 


TO  GET 
DIVIDE 


angstrom 

atmosphere  (normal) 

bar 

barn 

British  thermal  unit 
(thermochemical) 
calorie  (thermochemical) 
cal  (thermocheaicaD/cm2 
curie 

degree  (angle) 
degree  Fahrenheit 
electron  volt 
erg 

erg/second 

foot 

foot-pound-force 
gallon  (U.S.  liquid) 
inch 
Jerk 

Joule/kilogram  (J/kg) 

(radiation  dose  absorbed) 
kilo cons 
kip  (1000  lbf) 
kip /inch2  (ksl) 
ktap 

micron 

mil 

mile  (international) 
ounce 

pound-force  (lbf  avoirdupois) 
pound-force  inch 
pound-force/inch 
pound-force/foot2 
pound-force/inch2  (psi) 
pound-mass  (lbm  avoirdupois) 
pound-mas s - f oo c 2 
(moment  of  Inertia) 
pound-mass /foot2 

rad  (radiation  dose  absorbed) 
roentgen 


1.000  000  X  E  -10 
1.013  25  X  E  *2 
1.000  000  X  E  *2 
1.000  000  X  E  -28 
1.054  350  X  E  >3 


4.184 

4.184 

3.700 

1.745 

V(t* 

1.602 

1.000 

1.000 

3.048 

1.355 

3.785 

2.540 

1.000 

1.000 


000 

000  X  E  -2 
000  XE*1 
329  X  E  -2 
f*459.67)/1.8 
19  X  E  -19 
000  X  E  -7 
000  X  E  -7 
000  X  E  -1 
818 

412  X  E  -3 
000  X  E  -2 
000  X  E  *9 


shake 

slug 

torr  (a 


Hg,  0*C) 


4.183 

4.448  222  X  E  +3 
6.894  757  X  E  *3 
1.000  000  X  E  *2 

1.000  000  X  E  -6 
2.540  000  X  E  -5 
1.609  344  X  E  *3 
2.834  952  X  E  -2 
4.448  222 
1.129  848  X  E  -1 
1.751  268 
4.788  026  X  E  -2 
6.894  757 
4.535  924  X  E  -1 
4.214  011  X  E  -2 

1.601  846  X  E  +1 

1.000  000  X  E  -2 
2.579  760  X  E  -4 

1.000  000  X  E  -8 
1.459  390  X  E  *1 
1.333  22  X  E  -l 


meters  (m) 

kilo  pascal  (kFa) 
kilo  pascal  (kPa) 
meter2  (m2) 
joule  (J) 

joule  (J) 

mega  joule/m2  (MJ/m2) 
giga  becquerel  (GBq)* 
radian  (rad) 
degree  kelvin  (K) 
Joule  (J) 

Joule  (J) 
vatt  (W) 
meter  (a) 

Joule  (J) 
meter2  (m2) 
meter  (m) 
joule  (J) 

Gray  (Gy)** 

ter a joules 
nevton  (N) 
kilo  pascal  (kFa) 
nevton-second/m2 
(N-s/m2) 
meter  (m) 
meter  (a) 
meter  (m) 
kilogram  (kg) 
nevton  (N) 
nevcon-meter  (N*m) 
nevcon/meter  (N/m) 
kilo  pascal  (kPa) 
kilo  pascal  (kFa) 
kilogram  (kg) 
kilogram-meter2 
(kg*m2) 

kilogram/mecer2 

(kg/m2) 

Gray  (Gy)** 
coulomb /kilogram 
(C/kg) 
second  (s) 
kilogram  (kg) 
kilo  pascal  (kPa) 


*  The  becquerel  (Bq)  is  the  SI  unit  of  radioactivity;  1  Bq  ■  l  event/*. 
**The  Gray  (Gy)  is  Che  SI  unit  of  absorbed  radiation. 
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SECTION  1 


CAPABILITIES  OF  FRAM 


FRAM  (  No  acronym)  is  a  computer  code  which  computes  the  pressure  and  the  velocity 
of  a  particle  in  a  fluid  half  space  where  cavitation  may  occur.  The  cavitation  is  assumed 
to  result  from  an  exponentially  decaying  compression  wave  which  travels  to  the  surface 
and  is  reflected  as  a  tension  wave.  Because  the  water  does  not  sustain  tension  a  cavitated 
region  is  formed.  The  cavity  subsequently  closes  under  gravity  (lower  closing  shock)  and 
accretion  of  vapor  particles  (upper  closing  shock). 

The  model  of  the  fluid  and  the  mathematical  method  used  to  determine  the  boundaries 
of  the  cavitated  region  are  described  in  ref  [lj.  FRAM  computes  the  fret- field  solution,  that 
is,  the  field  variables  when  no  structure  is  present.  The  problem  of  the  interaction  with  a 
floating  structure  requires  the  use  of  an  approximate  interactive  scheme  such  as  the  plane 
wave  approximation  (P.W.A.)  or  the  doubly  assymptotic  approximation^.  A.  A.)  in  ref.  [l]. 

FRAM  does  not  treat  bottom  reflections  and  requires  the  explosion  to  be  located  a 
great  distance  away  from  the  observation  point  in  order  to  approximate  the  spherical 
waves  by  plane  waves.  It  treats  the  incident  wave  as  plane  and  steady  state  with  respect 
to  a  system  of  coordinates  moving  with  the  wavefront. 

The  one  dimensional  characteristics  method  is  used.  It  is  applied  to  the  two  di¬ 
mensional  problem  by  transforming  coordinates  ( X,Y,Z,T )  into  (y,  f)  where  y  *  Z,(  • 
y/lx7  +  Y7)  +  sT.  Combining  two  or  three  parameters  into  £  is  done  under  the  steady- 
state  assumption  that  the  cavity  is  unchanged  with  respect  to  a  coordinate  system  moving 
with  the  apparent  wave  speed  («).  Conservation  of  momentum  and  continuity  of  velocity 
across  the  closing  shocks  together  with  characteristic  relations  form  the  governing  system 
of  equations  . 

The  analytic  solution  provided  by  FRAM  has  the  advantage  over  numerical  methods 
that  is  that  it  requires  comparatively  small  computer  resources  and  sharper  wave  fronts  are 
obtained.  The  shortcomings  result  mainly  from  the  inability  to  include  bottom  reflection 
effects  and  the  divergence  and  decay  associated  with  propagation  away  from  the  explosive 
source. 
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SECTION  2 


DESCRIPTION  OF  THE  CODE 

2.1  PROGRAM  ORGANIZATION. 


Given  the  parameters  of  the  explosion  and  of  the  fluid,  the  first  step  consists  in  finding 
the  geometry  of  the  cavitated  region.  Along  with  the  geometry,  pressure  and  velocity  on 
the  boundary  are  also  computed.  The  second  step  takes  the  coordinates  (y,  f)  of  a  point  of 
observation  in  the  fluid  and  computes  its  pressure  and  velocity.  This  is  repeated  for  each 
time  increment  to  construct  a  time  history. 

The  main  program  calls  the  subroutine  CAVITY  which  computes  the  geometry  of  the 
cavitated  region.  It  is  formed  by  three  curves:  AB,  BC,  AC  (see  figure  1,  section  3  ). 
Having  charted  the  geometry  and  physical  properties  of  the  free-field,  one  can  proceed 
with  the  determination  of  the  properties  of  a  point  of  interest.  Subroutine  ECHO  takes 
the  coordinates  of  the  point,  finds  its  location  with  respect  to  the  cavity  using  subroutine 
ZONE,  and  finally  computes  pressure  and  velocity  at  that  point.  Finally  functions  G  and 
GG  correspond  to  the  characteristic  functions  g  and  to  its  derivative  /  along  the  Y  axis. 

There  are  two  types  of  parameters  that  must  be  specified:  those  pertaining  to  the 
material  properties  (e  speed  of  sound  in  water  ,  Pa  atmospheric  pressure,  to  fluid  density,  g 
gravitational  accelerarion)  and  to  the  explosion  (L  decay  length, ps  peak  pressure,  p  decay 
constant)  ,  and  the  coordinates  (y,  $)  of  the  points  of  observation.  All  parameters  and 
variables  are  given  in  a  non-dimensional  formulation.  In  the  example,  U.S.  customary 
units  are  used.  The  output  of  the  program  consists  in  the  curves  AB,AC,BC  forming  the 
cavity,  and  in  time  histories  of  p,  m,  n. 
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2.2  DESCRIPTION  OF  SUBROUTINES. 


The  present  version  of  the  code  is  designed  to  interact  with  the  finite  element  code 
SAP.  A  model  for  &  ship  provides  the  location  of  the  wet  nodes  at  which  pressure  and 
velocity  histories  are  computed  and  input  into  SAP.  The  SAP  input  is  read  from  tape  14 
and  written  on  tape  9.  It  is  not  difficult  however  to  extract  those  subroutines  relevant  to 
the  free-field  only  for  use  of  FRAM  with  other  codes.  The  subroutines  interacting  with 
SAP  are  READ  IT,  EMPER,  RESTWR. 

2.2.1  FRAM. 

The  main  program  calls  READIT  where  some  of  the  parameters  are  defined  and  the 
coordinates  of  the  point  of  observation  are  provided  by  SAP.  NWET  is  the  number  of 
wet  nodes,  KTIM  is  the  number  of  time  history  points.  The  explosion’s  parameters  are 
obtained  in  part  from  empirical  data  in  EMPIR.  CAVITY  is  called  for  each  point  of  obser¬ 
vation  and  ECHO  is  called  for  each  time  step.  The  output  is  written  for  SAP  in  RESTWR. 

Although  running  the  code  does  not  require  a  complete  understanding  of  the  theory, 
references  to  equations  and  figures  in  ref.[l]  have  been  included  for  the  sake  of  completeness. 
The  correspondence  between  some  of  the  symbols  in  ref.  [1]  and  in  FRAM  is  as  follows: 

CC=e,  PS=  p5,  L=L,  PA=  pA,  GAM=  io,  GV=  g  ,  T=  0,  TT=  fi,  AA=  a,  AL=  a l, 
YB=te,  SB=£b,  YC=y c,  SC=£c,  YA=yA,  SA=U,  FH=£,  F=%t,  XI(l)=Xfl,  XI(2)=Xa, 
R=R,E=A£,  A— A,  AS— fa,  GR=yr,  ST— £,  HS=A£,  YX=y,  PRE=/>,  VEL=m,  UEL=n, 
F -F,  G=G. 

2.2.2  CAVITY. 

In  the  following,  references  to  equations  and  figures  found  in  ref.[l]  will  appear  in 
parentheses.  This  subroutine  starts  by  defining  some  parameters  in  order  to  use  a  non- 
dimensional  formulation.  The  coordinates  of  point  B  are  given  by  finding  the  first  point 
on  the  opening  interface  AB  where  the  density  has  a  change  in  sign  due  to  closing  of  the 
cavity  (eq  (39)).  The  coordinates  of  point  A  are  found  on  AB  where  the  pressure  first 
drops  to  zero  (eq(51))  using  subroutine  NEWTON.  The  slopes  of  AC  and  BC  at  points  A 
and  B  are  also  evaluated  usingTaylor  series  expansions  (eqs(90)  and(83)). 

Next,  curves  AB  and  BC  are  constructed  incrementally,  with  an  increment  size  E 
which  is  determined  such  that  £s  coincides  with  the  £  value  of  one  of  the  discrete  points  of 
curve  AC.  The  curves  AC  and  BC  are  constructed  one  point  at  a  time.  Starting  from  the 
coordinates  of  a  point  I  and  the  slope  at  that  point,  one  can  determine  the  coordinates  of  the 
next  point,  1+1  and  the  slope  of  the  next  increment  by  a  first  order  linear  approximation. 

The  array  D(I,J)  contains  the  properties  of  the  curve’s  discrete  points.  J=1  corresponds 
to  curve  BC  or  lower  curve  and  J=2  to  curve  AC  or  upper  curve.  I  is  the  index  of  the 
discrete  points;  there  are  KK  more  points  on  AC  than  on  BC.  D(I,1)  is  the  characteristic 
//  (eq(97)),  D(I,2)  is  the  characteristic  gt  (eq(99)). 


Additional  arrays  which  do  not  appear  in  ref.  [1]  are  M(I)  and  P(I).  These  represent 
M=-y  +  P=v  +  6(  and  are  used  to  determine  the  location  of  the  point  of  observation 
with  respect  to  the  cavity,  i.e.  between  which  of  the  discrete  points  of  the  cavity  the 
characteristic  functions  defining  the  properties  of  the  point  of  interest  fall. 

S  is  augmented  at  each  step  by  an  increment  and  Y(I)  is  the  result  of  the  previous 
cycle  of  computations(eq(69)).  For  the  upper  curve,  the  characteristic  emanating  from  I 
is  reflected  on  the  free  surface  at  U  and  intersects  the  y  axis  at  T  or  the  curve  AC  itself 
at  R  see  (fig.(8)in  refjl]).  For  the  lower  curve,  R  is  defined  as  the  intersection  of  the 
characteristic  emanating  from  point  I  and  the  y  axis  (fig.(7.a)  in  ref[l]  ). 

The  upper  curve  is  constructed  first,  it  is  more  complex  than  the  lower  one  since  the 
characteristics  are  reflected  on  the  free  surface  and  fall  back  on  the  portion  of  the  curve 
freshly  constructed  at  point  R.  The  properties  at  point  1+1  are  deduced  from  those  at  point 
I  using  conservation  equations  and  geometry  (eqs(69)  and  (71))  in  which  the  expressions 
for  fu.tw.gw  are  needed.  If  R  falls  on  the  y  axis, then  it  corresponds  to  point  T  in  fig.(7.c) 
of  ref.  [l],  fu  is  determined  from  geomtricalrealtions  (eq(73)).  If  R  falls  on  the  beginning 
of  AC  ,  then  an  interpolation  between  the  two  adjacent  points  is  required  to  evaluate  fu. 
The  remaining  two  terms,  gw,g,w  are  simply  found  by  calling  the  functions  g with  the 
argument  yw-  g  is  given  by  its  definition  (  eqs(27),(33),  and(34)  for  a  point  on  the  y  axis). 

The  construction  of  the  lower  curve  requires  the  evaluation  of  a  similar  equation 
expressing  the  conservation  laws  and  geometrical  relations  (eq(68)),  and  in  particular 
gw,g’w,gr  which  is  simply  done  by  calling  the  functions. 

This  procedure  is  repeated  until  the  two  curves  intersect  at  point  C. 

2.2.3  ECHO. 

In  order  to  obtain  p,m  at  y,(  it  is  sufficient  to  know  f,g  (eqs.(102),  (103)).  Using 
geometry  and  characteristic  relations,  the  characteristics  at  X  are  expressed  in  terms  of 
similar  characteristics  at  points  where  these  are  known.  Fig. (8)  in  reffl]  shows  the  possible 
zones  in  which  the  point  may  fall.  A  zone  is  the  ensemble  of  points  with  same  expressions 
for  the  generalized  Riemann  invariants  (Table.  1).  ECHO  calls  ZONE  where  the  zone 
number  is  found  and  pressure  and  velocity  axe  computed  accordingly.  Results  are  converted 
in  dimensional  quantities  and  stored  in  0(NT,J)  for  plotting. 

2.2.4  ZONE. 

Z  is  the  zone  number,  F  and  G  are  defined  in  Table  1..  The  coordinates  of  the  points 
defined  in  Table  2.  are  computed  first.  Then  a  test  is  performed  to  determine  whether 
X  belongs  to  a  given  zone.  F  and  G  are  assigned  appropriate  values  which  allows  the 
computation  of  p,  m.  In  some  cases  it  is  necessary  to  call  SEARCH  to  locate  X. 

2.2.5  SEARCH. 

Using  the  arrays  P,  and  M  created  in  CAVITY  this  routine  finds  between  which  points 
of  AC  or  BC  the  characteristics  emanating  from  X  falls.  HFG  is  either  F  or  G.  When  X 


SECTION  3 


EXAMPLE 

Examples  of  output  are  presented  in  the  following  figures.  An  example  of  the  cavitated 
region  is  shown  in  Fig.l. 

The  location  of  the  charge  relative  to  the  submerged  surface  of  the  ship  and  origin  of 
SAP  coordinates  appear  in  Fig.2.  The  free-field  values  for  the  submerged  nodes  8,  18,  and 
28  which  are  specified  within  SAP  as  nodes  33,  43,  and  53  appear  in  Fig.3. 

Fig. 3  shows  in  (a)  the  free-field  pressure  which  is  characterized  by  a  pressure  pulse 
followed  by  a  surface  relief  effect  (pressure  drop  below  zero  line)  and  terminated  by  a 
secondary  pulse  due  to  the  closure  of  the  cavitated  region.  The  integral  of  the  pressure  is 
the  impulse  which  is  plotted  in  (b).  One  can  follow  in  (c)  the  correlation  between  vertical 
free-field  velocities  and  pressures  of  (a).  The  initial  upward  velocity  is  due  to  the  first 
spike  in  (a) .  A  free  fall  under  gravity  inside  the  cavitated  region  follows  and  is  terminated 
due  to  the  increase  in  pressure  resulting  from  the  cavitation  closure  (second  spike  in  (a)). 
These  free-field  velocities  are  integrated  to  give  displacements  shown  in  (d) . 
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